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■ We have calculated the complete one-loop corrections to the sfermion pair pro- 



duction process e"'"e~ — > fifj (/ = t,h,T,Vr) in the Minimal Supersymmetric Stan- 
dard Model. Our results also include the previously calculated SUSY-QCD cor- 
rections. We present the details of the renormalization scheme used. It is found 
that the weak corrections are of the same magnitude as the SUSY-QCD correc- 
tions at higher energies (y^ ~ ITeV). At these energies an important part of the 
weak corrections stems from the box contribution. This is best seen in sneutrino 
production. 



1 Introduction 



If supersymmetry (SUSY) is realized in Nature there should be two scalar particles 
(sfermions) fi, Jr corresponding to the two chirality states of each fermion /. The 
sfermions of the third generation play a special role as fi and may strongly mix 
(proportionally to the fermion mass), forming the two mass eigenstates /i and /2 (with 
f = t,b,T) . As a consequence one eigenstate (/i) can have a much lower mass than the 
other one. 

Sfermion pair production in e~^e~ colhsions, e~^e~ — > fi fj, = 1,2), has been studied 
extensively phenomenologically [J. The strong interest in sfermion production is mainly 
due to the fact that it gives access to one of the fundamental SUSY breaking parameters 
Af, the trilinear coupling parameter. It is clear that in e~^e~ —>■ ti tj and 6j bj gluon 
radiation and gluon exchange play an important role [21 E]- The SUSY-QCD corrections 
to these processes due to gluino and squark exchange were calculated in jH E| and found 
to become effective at ^/s > 500 GeV. Yukawa corrections [Hj were shown to be non 
negligible either. Very recently, while we were already working on the calculation of the 
full one-loop corrections to sfermion pair production within the Minimal Supersymmetric 
Standard Model (MSSM), such a calculation was presented in jZj. In this context, it is 
worthwhile to mention that the complete one-loop corrections to selectron and smuon pair 
production from threshold to high energies were calculated in jH]. 

For the calculation of higher order corrections, renormalization of the MSSM with an 
appropriate fixing of the SUSY parameters is necessary. Essentially, two methods were 
proposed in the on-shell scheme, one in [HI E]] and the other one in |1H IT^ for a review 
see [13] . Of course, both should lead to the same results for observables as masses, cross- 
sections, widths , etc.. The schemes differ in the fixing of the counterterms of some of 
the SUSY parameters as Mi,M2,/x, etc. Therefore the meaning of these parameters is 
different at loop level. However, at one-loop, in sfermion pair production this difference 
only matters in selectron pair production, e"^e~ — > Cj , and sneutrino pair produc- 
tion, e~^e~ ^ VeVe 1 due to the neutralino or chargino exchange being already present 
at tree level. Here we fix the SUSY parameters entering the sfermion mass matrices in 
the corresponding sfermion sector so that one has not to take into account any shifts in 
these parameters at one-loop level. In our case, we also have a different fixing of the fine 
structure constant a taking a{mz) as input, in contrast to [Zj where the Thomson limit 
is used according to [Hj. 

In this paper, we also present a full one-loop calculation within the MSSM for e"'"e~ — > 
fifjyf = t,b,T,h'T-. We compare our results with those obtained in [7]. Due to the 
complexity of such a calculation, an independent computation seems appropriate. We 
have calculated all graphs analytically and have written our own computer program for 
the numerical evaluation. Checks have been performed using the computational package 
fSlE]- In addition to the comparison with [7|, we studied different physical scenarios 
in our numerical analysis. We also include a study of the tau-sneutrino production which 
was not presented elsewhere. 
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2 Tree level 



The sfermion mixing is described by the sfermion mass matrix in the left-right basis 
the mass basis (/i, /2), f = t,b or t [T7| ITH]. 

afmf\ , i/m? 



^ \afmf mj^ J ^ ^ \ mf^ J 
where is a 2 x 2 rotation matrix with rotation angle 9j, which relates the mass 



eigenstates /i, i = 1,2, {rrij^ < rnj^) to the weak eigenstates fa, a = L,R, by fi = Rfafa 
and 



= M^^Q + {if-ej sm^e^) cos2/5mJ + mj , (2) 

m?^ = M^j^ ^1 + e/sin^6'vKCOs2/3m| + , (3) 

a/ = A/- -/i(tan/5)~^^''^ . (4) 

Mq, Ml, Mfj, Mf) and are soft SUSY breaking masses, ^/ is the trilinear scalar 
coupling parameter, fi the higgsino mass parameter, tan/3 = ^ is the ratio of the vac- 
uum expectation values of the two neutral Higgs doublet states , denotes the third 
component of the weak isospin of the fermion /, ej the electric charge in terms of the 
elementary charge e, and 6^ is the Weinberg angle. 

The mass eigenvalues and the mixing angle in terms of the primary parameters are 



m| = - (m| + m| =F —m% Y + 4a^m^ ) (5) 

/1, 2 2 ^ JR \ ^ fL fn/ J J J ^ ' 



^°^^/ = J ^/ , (0<^/<vr), (6) 



and the mass of the sneutrino is given by m^^ = M| + | m| cos 2(3. 
The tree-level cross-section of e~^e~ fifj is given by 



16 TT 



where 



-'77 ~ ^2 2^ ^ ■f^'' ' 



Asis — 



T-yZ = --^-^-^{ClKl + CrKr), (9) 
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and z) = {x — y — zY — ^yz. 

Here we use K^^i and C^^/j as the left- and right-handed couphngs of the electron to the 
photon and Z boson, respectively, 

Kl = Kr = 1, Ci = -^ + 4, CR = slr. (11) 

The matrix elements a|j come from the coupling of Zfifj, 

J _ ( 4(Jf cos^ ej - sl^Cf) -21 f sin 29 j \ 
""'^ ~ \ -2lf sin 2^; 4(/f cos^ - sl^ej) J' 

3 Radiative corrections 

The one-loop (renormalized) cross-section a"^™ is expressed as 

^■•^"(e+e- ^ fjj) = a'''' + Aa'^^^ + Aa^^, (13) 
where the symbol A denotes UV-finite quantities. 

As mentioned in the Introduction, the SUSY-QCD corrections (Ao"'^'"^) has already been 
calculated. In this paper, we give the result for the complete one-loop electroweak correc- 
tions (Acr^^) using the SUSY invariant dimensional reduction (DR) as the regularization 
scheme. The calculation was performed in the 't Hooft-Feynman gauge, ^ = 1. 
The electroweak corrections can be split further into the following UV-finite parts as: 

where Ao"^® and Aa^^ stand for the renormalized electron and sfermion vertex, A(jP™p 
for renormalized propagators and Aa^°^ for the box contribution. 
The renormalized electron vertex has the form 

/V^ K^(s, m| , m| ) , . 
Aa^^' = ^^-^-^(AT|,^ + AT|/ + ATJ|), (15) 



where 



ATI/ = '^^i^iAe,K, + AenKn), (16) 



2 2 / j: 

AT|/ = - ^f ""^y^ {AclCl + AcrCr + Aa^KL + Aa^K^), (17) 

ATJI = i6(s-m|)^ {Aa,C, + AanCn). (18) 

^^L,R and Aa^^ji consist of 3 parts, 

^eL,R = 5e% + 5et}, + 6e%, (19) 

AaL,R = 6a% + 5atl + 5a^l]^. (20) 



^^L^Ry ^^L^R correspond to the vertex corrections in Fig. [H 5ei\, ^cl^lr the wave- 
function corrections (Fig. ^ and Se^l^^, ^cl'lr correspond to the counterterms. 
The renormahzed sfermion vertex has a similar form, 

Acr^f = ^ (at;/ + AT// + ATj/) , (21) 



16 TT 



where 



AT// = '^^^^^{Kl + KD, 



77 g2 



(22) 



AV = - , /^ 2 , ^KlCl + i^«g^) ((Ae;),,a4. + ^.,(Aa^).,) , (23) 

= '^^^^^.{Cl + Cl). (24) 
16 (s — m^j 

{Aef)ij and {Aaf)ij can also be split into vertex corrections (Fig. wave-function cor- 
rections and counterterms, 

(Ae/).. = (5e;)S;) + (5e;)S/ + (5e;)g) , (25) 
(Aa/)., = {Saf)\f + {6afp + {6af)\f . (26) 

Explicit formulae for the remaining contributions Aa^^"'^ and Aa^°^ along with the for- 
mulae for g^l'r'^^ (^L^R'^\^^f)i]'^'^\ i^(^f)ij'^''^^ will t)6 given elsewhere. 

3.1 Fixing of the parameters 

The parameters occurring at the tree level, i.e. m^y, m^, sin 6'v^- are fixed as in ^Hl and 6j 
is fixed according to 120] • The difference to [7] is our fixing of the electric charge since 
we use as input parameter for a the MS value at the Z-pole, a = a{mz)\jj^ = e^/(47r). 
The counterterm then is ^Tl 



6e 1 e 



e (47r)2 6 



4EiV^ej(A + log^)+EEiV.^4(A + log-^ 
+ 4 ^ f A + log ^) + f f A + log - 22 f A + log 



with = m^j \/ rrif < rnz and Xt = rrit. N[, is the colour factor, A'^^ = 1,3 for (s)leptons 
and (s)quarks, respectively. A denotes the UV divergence factor, A = 2/e — 7 -|- log47r. 

3.2 Real corrections 

The cross-section a{e~^e^ fifj) is IR-divergent owing to the photon mass being zero. 
This is remedied by introducing a small mass A and including also the Bremsstrahlung 
i.e. a(e+e-^/47) (Fig. ED. 
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+ crossed graph 




+ crossed graph + crossed graph 




+ crossed graph 



Figure 1: Vertex and box diagrams relevant to the calculation of the electroweak correc- 
tions to e+e" —>■ fify 
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Figure 3: Bremsstrahlung diagrams relevant to the calculation of the real corrections to 
Summing those two contributions yields an IR-finite result for the physical value A = 0, 

In our calculation, we used a soft-photon approximation ^3] which reproduces the diver- 
gence pattern correctly but introduces a cut AE on the energy of the radiated photon. 
The explicit formulae will be given elsewhere. 

4 Outline of the calculation 

The results presented in this paper come from a full analytic calculation where we have 
neglected the electron mass (mg = 0) except for the QED corrections. 
In the case of this particular process, we can separate off the QED corrections on the 
basis of Feynman diagrams in a gauge invariant way. The QED corrections consist of all 
the diagrams that have an additional photon added to the tree level and therefore also 
include the whole real corrections. The weak corrections are then UV and IR finite and 
AE independent. All numerical results show only the weak corrections as the QED part 
is very sensitive to AE. 

The numerical calculation was performed using the packages LoopTools and FF ^Hl- The 
results were checked against the results for the Yukawa approximation presented in jB] 
where our results match except for a minus sign already pointed out in [7] . 
Furthermore, we did our own independent calculation based on FeynArts and FormCalc 
packages ^Hl checking all individual Feynman diagrams at the amplitude level. In addi- 
tion, we used the packages to check the complete results using the same renormalization 
scheme as in the analytical calculation. 

We also compared our results with [7j where we were able to reproduce all the results 
apart from minor differences due to the different renormalization of the fine structure 
constant. 

5 Numerical analysis 

In the following numerical examples, we assume Mq = Mq^ = yiVfj^^ = ^Mfj^ = M^^ = 

= Mn = Mfr = Mn = Mf = Mfp for the first, second and third gener- 
ics VI, 2 (Jl,2 i^l,2 1^1,2 i^'1,2 ' ° 

ation soft SUSY breaking masses, and all trilinear couplings are set to a common value 
A. For the standard model parameters we take ruz = 91.1875 GeV, rriw = 80.45 GeV, 
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sin^^^^ = 1 - m^/m|, a = 1/127.934, = 174.3 GeV, and rrib = 4.7 GeV. M' is fixed 
M by TTig = {as{'mg)/a) sin^ 9wM. 

Below we show plots for three different scenarios. On the left, there are the total and 
tree-level cross-sections for all channels e~^e~ —>■ fifj, f = t,b,T, i = 1, 2, and on the right 
we picked out one of the channels and show a separation into the convergent subclasses 
described in the text above. 

Care has to be taken when comparing plots for different sfermions in the same scenario 
because we always fix the SUSY parameters (Mq,/^ etc.) in the sector of the produced 
sfermions. 

One can see that the total corrections in squark production are dominated by SUSY-QCD 
where the biggest contribution comes from the gluon. The gluino part is small compared 
to the gluon one as already shown in jS]. At higher energy (1 TeV — 1.5 TeV), other 
corrections can grow to a size comparable to the SUSY-QCD contribution. The leading 
weak contribution at high energy comes from the box diagrams. As can be seen in Fig. 
mini the box diagrams give a negative contribution rising with energy. This can be checked 
against the Sudakov approximation at high energies [221 121] ? where the box contributions 
give the leading correction. The sum of the other two corrections, the vertex and the 
propagator corrections, is small compared to the box contribution. 

The weak corrections computed here have a significant effect in the high energy region, 
and in the case of squarks act against the SUSY-QCD corrections. For stau production 
they can go up to -10% in certain cases. 

In addition, we also show the tau-sneutrino production. For detection, it is necessary that 
not only the decay channel to Xi is open. Our three scenarios allow decays into charginos. 
Sneutrino production is the only case where not only the SUSY-QCD corrections are not 
present but also the Yukawa corrections are small. This is due to the fact that the dia- 
grams including a neutralino or a neutral Higgs boson in the loop are missing (in Yukawa 
approximation). Therefore, the sneutrino production shows a particular dominance of the 
box corrections. In Fig. [7|one sees that the tree level is almost identical in two of our 
scenarios due to the small difference in the sneutrino mass. The total cross-section in the 
two scenarios is also very similar as the different vertex corrections and the propagator 
corrections are together below 5% and the boxes give the leading contribution. 

6 Conclusion 

We have calculated the complete one-loop corrections to stop, sbottom, stau and tau- 
sneutrino production. The calculation was performed in an analytical way with an in- 
dependent check using the FeynArts and FormCalc computer packages ^H]- Our way of 
fixing the fine structure constant a gives a higher tree-level cross-section and therefore 
smaller radiative corrections compared to [7j. The corrections are typically of 5-10% and 
thus not negligible at a future linear collider. 
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SCENARIO 1- gaugino: 



The parameters are set to 

{M, /i, A, tan /3, m^, Mq] = {200 GeV, 1000 GeV, -500 GeV, 20, 300 GeV, 400 GeV} 
The masses of the sfermions in this scenario are 

nii^^^ = {276, 520} GeV m^,^ ^ = {285, 526} GeV m^, ^ = {354, 446} GeV 




600 800 1000 1200 1400 1000 1200 1400 



^ [GeV] ^ [GeV] 



Figure 4: Scenario 1: (a,b,c) Total (sohd) and tree- level (dashed) cross-sections as a 
function of ^/s ; (d,e,f) Relative corrections to one of the channels from plots (a,b,c) 
split into convergent subclasses, (thick solid- Aa^°^^^, solid- Aa^*"^, dashed- Aa^°^, 
dash-dotted- Aa™''*^^, dash-dot-dotted- ActP^p) 
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SCENARIO 2- higgsino: 



The parameters are set to 

{M, /i, A, tan /3, m^, M^} = {800 GeV, -200 GeV, 250 GeV, 30, 300 GeV, 250 GeV} 
The masses of the sfermions in this scenario are 

{200, 360} GeV m^^ ^ = {203, 318} GeV ruf^^^ = {231, 275} GeV 
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Figure 5: Scenario 2: (a,b,c) Total (sohd) and tree-level (dashed) cross-sections as a 
function of ^/s ; (d,e,f) Relative corrections to one of the channels from plots (a,b,c) 
split into convergent subclasses, (thick solid- Aa^°^^^, solid- Aa^*"^, dashed- Aa^°^, 
dash-dotted- Aa™''*^^, dash-dot-dotted- ActP^p) 
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SCENARIO 3- mixed: 



The parameters are set to 

{M, /X, A, tan /3, m^, Mq] = {200 GeV, 200 GeV, -800 GeV, 10, 300 GeV, 400 GeV} 
The masses of the sfermions in this scenario are 

{172,563} GeV = {398,446} GeV nif^^^ = {396,409} GeV 
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Figure 6: Scenario 3: (a,b,c) Total (sohd) and tree-level (dashed) cross-sections as a 
function of ^/s ; (d,e,f) Relative corrections to one of the channels from plots (a,b,c) 
split into convergent subclasses, (thick solid- Aa^°^^^, solid- Aa^*"^, dashed- Aa^°^, 
dash-dotted- Aa™''*^^, dash-dot-dotted- ActP^p) 
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SNEUTRINO PLOTS: 



Here we show the sneutrino production in all three scenarios. The masses of the sneutrino 
in the scenarios are 



394.795 GeV mj>^ = 242 GeV m^^ = 394.873 GeV 
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Figure 7: Sneutrino plots: (a,b,c) Total (solid) and tree-level (dashed) cross-sections 
as a function of ^/s ; (d,e,f) Relative corrections split into convergent subclasses, 
(thick solid- Aa*°*^', dashed- Aa^™, dash-dotted- Aa™''*''^, dash-dot-dotted- ActP^p) 
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